CO 

o 

CN 

Oh 
< 

in 

CN 



CM 

< 

c3 



> 



Boundedness and blow-up of solutions for a 
nonlinear elliptic system arising in probability 

and stochastic processes 

Dragos-Patru Covei 

Constantin Brancusi University of Tg-Jiu, Gorj, Romania 

E-mail: covei.dragos@yahoo.com 

Abstract 

The main objective in the present paper is to obtain the existence results for bounded and 
unbounded solutions of some quasilinear elliptic systems. Related results as obtained here have 
been established recently in [C. 0. Alves and A. R.F. de Holanda, Existence of blow-up solutions 
for a class of elliptic systems, Differential Integral Equations, Volume 26, Number 1/2 (2013), 
Pages 105-118.]. Also, we present some references to give the connection between these type 
of problems with probability and stochastic processes, hoping that these are interesting for the 
audience of analysts likely to read this paper. 
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^1 Introduction 

o ■ 

(ryfhe question of existence of solutions for elliptic equation of the form 

>■ \u = f (x, u) in Q, (1.1) 

• i— i ' 

rNwas studied by many researchers (see Bandle [2], the author [1], Lair [TT], Matero [2], Mohammed 
cd[TB] with their references). This work is devoted to the study of the more general nonlinear elliptic 




problems of the type 

F Ui (s,iti,...,iii,...,it d ) infi, 

(1.2) 

where d > 1 is integer, Q C M. N (N > 1) is a smooth, bounded domain or Q = M. N , A p Ui := 
div (|Vuj| p_ V«j) (1 < p < oo) is the p-Laplacian operator and F u . (i = l,...,d) stands for the 
derivatives of a continuously differentiable function F : Q x [R + ] d — > IR + in (u\, ...,Ud)- For the case 
Q = M. , we also consider the following class of elliptic systems: 

A p Ui = at (x) F Ui (x, Mi, ..., Ui, ..., u d ) in R N , 

^ > in R N , (1.3) 



where a« : WL N — > (0, oo) are suitable functions. Associated with the class of systems (II. 3p . our main 
result is concerned with the existence of entire large solutions, that is, solutions (ui, ■■■,u ( t) satisfying 
Ui (x) — > oo as |x| — > oo for all i = 1, ..., d. 

The interest on systems fll.2p - fll.3B comes from some problems studied in the works of Lasry-Lions 
[T2] , Busca-Sirakov [3] and Dynkin [5] where the authors give the connection between these type of 
problems with probability and stochastic processes and from the recently work of Alves and Holanda 
[1] where these systems are considered for the case d — p — 2, in terms of the pure mathematics. 
The difference between our work and the paper by [1] is that: our systems can have any number 
of equations, the potential functions a^ cover more general properties and that we use in the proofs 
theories for quasilinear operators instead of the theories for linear operators used by Alves- Holanda 

ID- 

To begin with our results we make the following convention: we say that a function h : [0, +oo) — > 
[0, oo) belongs to T if 

h G C 1 ([0,oo)), h(0) = 0, ti(t) >0WG [0,oo), 
h(t) > W G (0, oo) 

and the Keller- Osserman [10], [17] condition is satisfied, that is, 

1 



i H{ty n 



i//' 



dt < oo, 



where H (t) = f Q h (s) ds. 

Our main result for problem (II. 2p on a bounded domain is the following: 

Theorem 1.1. Suppose Q is a smooth, bounded domain in R N and that there exist f t , g G T 
satisfying 

F u (x, h, ..., U, ..., t d ) > fi (U) VieO, U > and i = l,...,d (1.4) 

and 

g (i) > max {F t (x, t, ..., t)} Vx G H, t > 0. (1.5) 

Then: 

1. problem hi. 2^1 admits a positive solution with boundary condition 



Ui = a>i on dil 
2 = 1, ..., d 
a.i G (0, oo) . 

2. problem U.2)) admits a positive solution with the boundary condition 



1.6) 



Ui = oo on dQ 



l,...,d 



;i-7) 



where Ui = oo on dfl should be understood as Ui (x) = oo as dis (x, dQ) — > 0. 

3. problem U.2)) admits a positive solution with boundary condition: there are io,jo £ {!)•••> d} 
such that 



u io = oo on dQ 



Uj < oo on dQ for any j ^ i 
and the set {1, ...,d} is crossed by i respectively j . 



;i.8) 



Our next result is related to the existence of a solution for system (11.31) . For expressing the next 
result, we assume that functions Oj (i = 1, ...,d) satisfy the following conditions: 

<H (x) > for all xeR N and a { G C%£ (R N ) , & G (0, 1) (1.9) 

and that the quasilinear system 

d 

—A p z (x) = Yl a i ( x ) f° r x e ^ JV > z ( x ) — >" as \ x \ — >" °° (1-10) 

4 = 1 

has a C 1 -upper solution, in the sense that 

J RN | Vzf" 2 Vz ■ W<pdx > J RN Z ti Oi ( x ) <f>dx, <P G C °° (R N ) , > 
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G C 1 (R N ) , 2 > in R N , z (x) ->■ as |z| ->■ oo. 



Theorem 1.2. Assume that (fl.^| )-( TOip /ioW. T/ien system U.3\) has an entire large C l -solution 
(in the distribution sense). 

To prepare for proving our theorems, we need some additional results. 

2 Preliminary results 

Let fi C R N (N > 2) be a smooth, bounded domain in R N and 1 < p < oo. The first auxiliary result 
can be seen in the paper of Matero [EJ pp. 233]. 

Lemma 2.1. Assume that g meets the conditions: g is a continuous, positive, increasing function 
on R + , and g(0)=0. Lethe W 1 * (fi) be such that (G o h) G L 1 (fi) w/iere G (s) = f° g (t) dt. Then 
there exists a unique u G W 1,p (fi) which (weakly) solves the problem 

A p u (x) = g{u (x)) , x G fi, 
u(x) — h(x) , x G <9fi. 

Furthermore, if u\ and u 2 are the solutions corresponding to hi and h 2 with hi < h 2 on dfl, then 
Ui < u 2 in fi. Finally, there exists an a G (0, 1) such that u G C 1,/3 (D) (f3 G (0, 1)) for any compact 
set D C fi. 

The following comparison principle is proved in the article of Sakaguchi [18] (or consult some 
ideas of the proof in work of Tolksdorf [T9], Lemma 3.1.]). 

Lemma 2.2. Let u, v G W l,p (fi) satisfy —A p u < —A p v for x G fi, in the weak sense. If u < v on 
dfl then u < v in Q. 

The following Lemma is essentially and can be found in Zuodong (|22j). 

Lemma 2.3. Suppose f G J 7 and that u G W l( £ (fi) fl C (fi) satisfies 

- f \Vu\ p ~ 2 VuV^>dx = f f(u) <pdx,V<p G C °° (fi) . 
Jn Jn 

Then, there exists a monotone decreasing function \i : (0, oo) — > (0, oo) determined by f such that 

u(x) < \x (dist (x, dfl)) \/x G fi. 

Moreover, 

limg (t) = oo, lim g (t) = — oo. 

t— >0 t— >oo 



Next, we begin with recalling the definition of sub and super-solution used in the present context. 
The system that we will study is the following 

A p Ui = G Ui (x, «i, ..., u h ..., u d ) in fi 

Ui = fi on <9fi (2-1) 

i — 1, ..., d 

where fi G W x,p (fi) and G (x, t±, ..., ti, ..., td) : fi x [R] -> R is measurable in x £ 11, continuously 
differentiable in ti G R, and satisfies the following condition: for each T, > fixed (z = 1, ..., d), there 
exists C = C (T^ > such that 

|G(x,* 1 ,...,i d )|<CV(x,ii,...,i d )Gfix [-T,,T,] d . (2.2) 

Now we introduce the concept of sub- and super- solution in the weak sense. 

Definition 2.1. By definition {u 1) ... ) u d ) G [W 1,p (12)] zs a (weak) sub-solution to Ii2. 1\) , if u { < fi 
on <9fi, 

/ |Vmj| p ~ VujVifidx — / G„. (x, u x , ■■■■ l Ud) 4>dx < 

/or a// G C^° (12) u>rt/i > and z = 1, ..., d. 

Similarly (ui,...,Ud) G [W 1,p (12)] zs a (weak) super-solution to \2. 1\) if in the above the reverse 
inequalities hold. 

The following result holds: 

Lemma 2.4. Suppose (u Xl ... 1 u d ) is a sub-solution while (ui,...,Tid) is a super-solution to problem 
Ii2. 1\) . and assume that there are constants g ±i ,a i G R such that 

On < Hi < ^i < Hi almost every where in 12. 

If \2. 2\) holds, then there exists a weak solution (ui,...,Ud) G \W 1,P (12)] of \2. 1\) . satisfying the 
condition 

u_i<Ui< Ui almost everywhere in 12. 

We will not give the proof here since he can now proved as in [TJ Theorem 2.1, pp. 110] with 
some ideas from 



3 Proof of main results 

In this section, we will prove the main results of this paper. 

3.1 Proof of Theorem [LJ] 
3.1.1 Proof of 1: 

In what follows, we denote by ip G W 1,p (12) the unique positive solution of the problem 

- / n | V?/f " 2 VipV(j)dx = J n g (tp) <j)dx in fi, V0 G C °° (fi) with > 0, 
ip > in fi 
ip = m on <9fi 



where m = min {a±, ..., a d }, which exists and minimizes the Euler-Lagrange functional 



on the set 



J(ij>)=: I l-\Vtp\ p + G{iP{x)))dx 



K = {v G L 1 (fi) \v-me W^' p (n) and (Gov)e L 1 (fi) } 



i.e., ip meets the boundary condition (ip — m) G W ' p (Q) in the weak sense [see Lemma [2.11 and [2J 
Paragraph 6, pp. 12]]. Then, 

— f n \^4>\ P ~ VipV(f)dx = J n g (ijj) <pdx > j n F^ (x,ip, — ,i>) 4>dx in Q, 
i\) = m < on on dQ 
i = 1, ..., d 

and so {u x , ..., u d ) — (i/;, ...,i/)) is a sub-solution. Clearly, (mi, ...,«d) = (M, ..., M) , with 

M = max{aj \i = 1, ..., d} , 



is a super-solution of fll.2)l with the boundary condition (II. 6p . We prove that, y^ < u l for all 

i = l, ..., d. Indeed, 

f 

— ApMj = — <? (mj) < — A p Ui = in O, 

Ui = m < Hi = M on dfl, 

i = 1, ..., d 

and then with the use of Lemma 12.21 it follows that y^ < w, in f2. According with the proof of the 
Lemma [2. 4| there exists a critical point (u±, ■ ■■,u d ) G \W 1,P (JT)] which minimize the Euler-Lagrange 
functional 



I(u x ,...,u d ) 



- / 2_.\^ u if dx + I F (x,u%, ...,Ud)dx 
P Jn i=l Jn 



that solve, in the weak sense, the system 



A p Ui = F Ui (x, Mi, ..., Ui, ..., u d ) in Q 
Ui = on on dfl 
i = 1 d 



(P a 



and satisfying ip < u^ < M in Q for all z = l,...,d. Since U{ G L^ c (Q), by the regularity theory 
(3 H31 ED], it follows that u.i G C l (O). 

3.1.2 Proof of 2: 

To study this case, we begin considering the system 



A p Ui = F Ui (x, ui, ..., Ui, ...,u d ) in Vt, 
Ui = n on dfl, 
i = 1, ..., d. 



(3-1] 



Then, by the finite case above, problem (13. ip has a solution (m™, ...,Uj). 

We prove that the sequence of solutions (w™, ...,u^) can be chosen satisfying the inequality 

< < < +1 for all i = 1, ..., d and n G N. (3.2) 

To prove this, we consider the solution (u\, ...,u^) of the problem 

A p Ui = F Ui (x, Mi, ..., Ui, ..., u d ) in Q, 

Ui = 1 on dft, (3.3) 

z = 1, ..., d. 
and note that it is a sub-solution of 

A p Ui = F Ui (x, Mi, ..., Ui, ..., M d ) in fi, 

^ = 2 on ^, (3.4) 

i = 1, ..., d 

while the pair (Mi, ..., Mi) is a super- solution for Mi = 2. Once < u.; (a;) < 2 (z = 1, ..., d) Wx G O 
, Lemma EH implies that there exists a solution (m^, ..., m^) of 

A^ = F Ui (x, ui, ..., Ui, ..., u d ) in fi, 
tij = 2 on <9fi, 

i = l, ..., d. 

satisfying u\ (x) < u\ (x). Using the argument above, for each M n — n+ 1; n — 1,2, ..., we get a 
solution (m™, ...,u^) °f (13 .11) . which is a sub-solution, and the pair (M n ,...,M n ) is a super-solution 
respectively of 

ApUi = F Ui (x, Ui, ..., Uj, ..., Ud) in fi, 

n,j = n + 1 on <9fi, , 

2 = 1, ..., <i. 

Thereby, the sequence of solutions (it™, ...,u2) satisfies the inequality (13. 2p . Finally, we construct an 
upper bound of the sequence. More exactly, we show that {(u™, ..., u^)} n>1 is uniformly bounded in 
any compact subset of fi. To this end, we begin recalling that by (11.41) 

ApWl > fi «) in fi 
u" > in fi 
u™ < n on <9fi 

with /j G J 7 . If uf (i = 1, ..., d) denote the unique solutions of the problems 

ApUi = fi (uj) in fi 
Ui > in fi 
Ui < n on <9fi 
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it follows from Lemma 12.21 that 

u? <u?mQ for all n > 1. 



By Lemma I2.3| there exist non- increasing continuous functions [i i : M + — > M + such that 

«F < Hi{dist(x,dQ)) Vn E N,Vx E fi and i = 1, ...,d 

showing that 

< «• (x) < < (x) < ^ (d (x)) VneN.ViGn (3.5) 

where d(x) = dist(x,dfi). Thus there exists a subsequence, still denoted again by u™, which con- 
verges to a function Mj in W 1,p (Q). In other words 

Ui (x) := lim w™ (x) for all x E f2 and i = 1, ..., d. 

The estimates ( 13. 5 p combined with the bootstrap argument yield that uf (x) — > Ui (x) in C 1 (K) 
for any compact subset K C CI. Furthermore, it is clear that, Ui (x) E C 1 (Q) and (u\, ■■■,Ud) is a 
solution of (II. 2p : that is, 

ApMj = F u . (x, Mi, ..., Mj, ..., Md) in Q, 
Mj > in Q, 
i = 1 d. 



., . .. , 



To complete the proof, it suffices to prove that (u\, ..., Ud) blows up at the boundary. Supposing for 
the sake of contradiction that Ui does not blow up at the boundary, there exist xo E dfi and (xfc) C CI 
such that 

lim Xk = Xo and lim Ui (xk) = Li E (0, oo) . 

k— )-oo fc— >oo 

In what follows, fix n > 4Lj and 5 > such that «" (x) > n/2 for all x E f^, where 

tt s = {x Ett \dist (x, an) < 6} . 

Then, for k large enough, Xk E fls and u" (x^,.) > 2Lj. Since 

< (a*) < < +1 (x fe ) < ... < < +J (a*) < ... < «< (x fc ) Vj, 

we have that U{ (xk) > 2 Li, which is a contradiction. Therefore, m« blows up at the boundary. This 
solution (ui,...,Ud) dominates all other solutions and is therefore commonly called blow-up/large 
solution. 

3.1.3 Proof of 3: 

Let (...,w™, ..., Uj o , ...) E C 1 (Q) be the solution of the problem (P a ) with a io = n, n E N, and a>j 
fixed. As in the previous case, the sequence u™ is bounded on a compact subset contained in Q, 
implying that there exist functions u io (i = l,...,d) satisfying u™ —> u io in C l (K) (i = l,...,d) for 
any compact subset K C Q. Moreover, the arguments used in the previous cases yield that Ui blows 
up at the boundary, that is, u™ = oo on dfl. Related to the sequence (uj ), we recall that 

A p w™ = F U . Q (x, <,..., m™, ..., ufj in fi, 
Uj = a jo on dil. 



Then, by the comparison principle w™ < aj Vx E fi and n > 1. Passing to the limit as n — > oo, we 
obtain that Uj < <x, for all igSl. 

Claim. Let x E <9f2 and (x fe ) C Q be a sequence with Xk — > x . Then w J0 (x fe ) — > Oj as k — ► oo. 



Indeed, if the limit does not hold, there exist e > and a subsequence of (xk), still denoted by 
itself, such that 

Xk —> Xq and Uj (xfc) < Oj — £ V/c G N. (3.6) 

Since Uj = Oj on <9f2 and is continuous, there is some 5 > such that w J0 (x fc ) > ct, — |, Vx G f^. 
Hence, for k large enough, x fc G Q$ an d n }„ > n j ( x fc) > «j — f > a io ~~ e wm °h contradicts (13. 6p . 
From this claim, we can continuously extend the function Uj from Q to Q by considering Uj (x) = aj 
on dQ, concluding this way the proof of the Finite and infinite case. 

3.2 Proof of Theorem IT21 



Firstly, we provide a sub-solution for the problem ( II. 3p . To do this we consider the function w 
M. N — > [0, oo) implicitly defined by 



z (x) 



-dt, x eW\ 



Note that weC 1 (R^, (0, oo)), w (x) — > +oo as |x| — > oo and 

Vz(x) = -g- 1/{p - 1] {w(x))Vw(x) 
|Vu> (x)\ p ~ Vw(x) = — g (w (x)) \Vz (x)\ p ~ Vz(x). 

Given G C£° (R^) , > we have 

|Vw(x)| p ~ Vw (x) V(j)dx = I — g (w (x)) |Vz (x)\ p ~ Vz{x)V<pdx 



(3.7) 
(3.8) 



/ div \g (w (x)) \Vz (x) \ p Vz (x)] 4>dx. 

JRN 



Computing the derivatives in the integrand of the expression just above, in the distribution sense, 
using (13. 7p we get, 



| Vu> (x) | p Vto (x) V(f)dx = I g (w (x)) A p z (x) 4>dx 

g' (w (x)) g~ (w (x)) \Vz (x) \ p <fidx. 
Using the fact that g G J 7 and that z (x) is an upper solution of (11.111) we derive the inequality 

/ |Vu> (x)\ p ~ Vw{x)V(j)dx — I g (w (x)) A p z (x) <pdx < 0, 

Jr n Jr n 

and so 

\Vw (x)\ p ~ Vw (x) V(f)dx + g (w (x)) > o, (x) I <f>dx < 0, 

Jn» \f^ J 

which together with (II. 5p leads to 

\Vw(x)\ p ~ Vw{x)V(j)dx> I Oj (x) F u . (x,w (x) , ...,w (x)) 4>dx 



for all i = 1, . .. , d. 

In the next, we consider the system 



A p Ui = di (x) F Ui (x, Mi, ..., u u ...,u d ) in B, n 



Ui = w n in dB n , 



i = l,...,d, 



(3.9) 



where B n is the open ball of radius n centered at the origin and w n = max l£ -g w(x). Clearly, 
(w,...,w) and (w n , ...,w n ) are a sub-solution and super-solution for (13.91) respectively. Thus, by 
Theorem 11.1] there is a solution (w™, ...,u2) G [W 1,p (B n )] of (13 .9p satisfying w (x) < uf < w n for 
all x G B n and i = 1, ..., d. By the comparison Lemma w™ < w™ +1 and moreover w™ G C 1 (-B n )- For 
m > 1 and n > m + 1 consider the family of systems 

A p < = <fi «) in B m+1 

i = 1, ..., <i, 

where a" = min xe ^ a^ (x) > 0. Arguing as in the previous sections, there are monotone decreasing 
functions u°; : (0, oo) — > (0, oo) determined by fi such that 

w (x) < u? (x) < nf {dist (x, dB m+x )) Vx € B m+l 

from which it follows that 

w (x) < m™ (x) < Mf for all n e N, x G 5 m , i = l,...,d 

and for some positive constants M™. As a consequence, there is Ui G L°° (B m ) (i = 1, ...,d) such 
that for some sub-sequence of uf, still denoted by itself, we get 

u" — > Ui (i = 1, ..., d) pointwisely in S m (Vm > 1). 

Therefore, («i, ..., Ud) G C 1 (IR^) and is a solution for the system 

A p Ui = at (x) F Ui (x, Wi, ..., Mj, ..., u d ) in M^ 

Ui > in R N 

% — 1, ..., <i 
v 

satisfying 

w (x) < «i (x) for all x G M and % = 1, ..., rf. (3.10) 

Letting |x| — > oo in (I3.10p it follows that (ui, ...,Md) is a large entire solution for (II. 3p . 

4 Remarks 

4.1 The probabilistic connection 

Assume that ip belongs to a wide class \I/ of monotone increasing convex functions. There is an area 
in probability theory where boundary-blow-up problem 

An = i/) (u) in fi 

u = oo on dfl 

arise (see the paper [0] or directly the book [7] for details). The area is known as the theory of 
super diffusions, a theory which provides a mathematical model of a random evolution of a cloud of 
particles. Indeed, given any bounded open set fl in the N-dimensional Euclidean space, and any 
finite measure u we may associate with these the exit measure from Q i.e. (Aq,P m ), a random 
measure which can be constructed by a passage to the limit from a particles system. Particles 
perform independently A-diffusions and they produce, at their death time, a random offspring (cf. 
[8]). Py, is a probability measure determined by the initial mass distribution u of the offspring and Xq, 
corresponds to the instantaneous mass distribution of the random evolution cloud. Then procedding 
in this way, one can obtain any function ip from a subclass \l/o of \& which contains u 1 with 1 < 7 < 2. 
Dynkin [B], also provided a simple probabilistic representation of the solution for the class of problems 
u 1 (1 < 7 < 2), in terms of the so-called exit measure of the associated superprocess. Moreover, the 
author say that a probabilistic interpretation is known only for 1 < 7 < 2. 



4.2 The stochastic connection 

In the next we present a simple example from the paper Lasry-Lions [T2] in order to connect our 
problems with stochastic processes. To do this, we consider the stochastic differential equation 

dX t = atdt + dB t 
X = x G Q 

where B t is the Brownian motion (on a standard probability space) and a t = a (X t ) is a feedback 
control. We also consider the class A of all (feedback) controls that keep the diffusion process X t 
inside the domain Q for any time t > a.s.. The criterion for optimality is given by the cost functional 
(E is the expected value, C q = (q — l)q~ q > and q' = -zr) : 

f(X t ) + C q \a(X t )f 1 e-Ht . 

assigned cost cost of the control I dlscount tactor 




Then 



inf J (x, a) 



where 

A = {a e C° (n,R N ) :X t eQ,Wt>0 a.s.} 

is achieved and defines the value function 

u (x) = inf J (x, a) 

that solves, in the case 1 < q < 2 (dynamic programming principle), the problem 

Au+ f (x) = u+ \Vu\ q 

(4.1) 

u (x) — > oo as d (x) — > 

where d (x) = dist (x, dQ) . 

We also remark from the work of Busca-Sirakov (3] that the above disscussion for the scalar 
equation (14. Xp suggest that the solutions of the system (1.1.) can be viewed as the value function of 
a stochastic control process, and the boundary conditions then means that the process is discouraged 
to leave the domain by setting an infinite cost on the boundary. For a more detailed discussion about 
practical applications where such problems appear we advise the reader the introduction of the work 

Hi- 
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